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ABSTRACT 

R-modes can generate strong magnetic fields in the core of accreting millisecond 
neutron stars (NSs). The diffusion of these fields outside the core causes the growth of 
the external magnetic field and thus it affects the evolution of the spin down rates P of 
the millisecond pulsars (MSPs). The diffusion of the internal magnetic field provides 
a new evolutionary path for the MSPs. This scenario could explain the large P of the 
pulsars J1823-3021A and J1824-2452A. 
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1 INTRODUCTION 

In recent papers it has been proposed a mechanism 
based on the r-mode instability for the evolution of 
the internal magnetic field in milli s econd accreting com- 
pact stars llCiiofano fc Dragol l2010l : iBonanno et al.l l201ll : 
ICuofano et al.l I2OI2I ). R-modes represent a class of oscil- 
lation modes in rotating neutron stars, which are unsta- 
ble with respect to the emi ssion of gravitational waves 
lAndersson fc Kokkotad l2001h . Only the existence of effi- 
cient damping mechanisms can suppress this instability and 
allows neutron stars to rotate at high frequencies. An im- 
portant mechanism for damping is dissipation which can be 
provided by shear and bulk viscosities converting the rota- 
tional energy of the star into heat. Another very efficient 
source of damping is given by the coupling of the r-modes 
with the pre-existing poloidal magnetic field of the star: it 
has been shown that r-modes are responsible for the forma- 
tion of huge internal toroidal magnetic fields which could 
reach values of 10^^ G or larger in the case of newly born 
neutron stars. Smaller values are obtained in the case of 
old and accreting neutron stars. A large fraction of the 
rotational e nergy of the star is thus stored in the ma g- 
netic field (iRezzqlla et all l200d: iRezzoUa et all I2OOI al ibi: 



ISa fc Tomel 120051 . I2OO6I : ICuofano fc Dragol I2OI0I ). This con- 
figuration becomes unstable due to the so-called Tayler in- 
stability (Braithwaitc 2009; Cuofano & Drago 201Q): a new 
poloidal component is generated which can then be wound- 
up itself closing the dynamo loop. The internal fields evolve 
into a stable configuration in which the toroida l component 
can b e significantly stronger than the poloidal (jBraithwaitd 
I2OO9D . 
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Here we analyze the evolution of the external magnetic 
field during/after accretion due to the diffusion outside the 
core of the magnetic field generated by r-modes. The in- 
crease of the external magnetic field implies a growth of the 
spin down rate P of the millisecond NSs. Finally we show 
that this mechanism can account for the large P observed 
for the two pulsars J1823-3021A and J1824-2452A. 



2 EVOLUTION OF THE EXTERNAL 
MAGNETIC FIELD 

The toroidal magnetic fields generated in the core of ac- 
creting millise cond neutron stars are in the range Bj"^ ~ 
[10^^-10"] G iGuofano fc Dragdl2010l : [Guofano et a"Lll2012l '). 
These strong toroidal fields can be stabilized with respect 
to the T ayler instabiUty by a much smaller poloidal com- 
ponent (jBraithwaita 120091 ') that we assume in the range 
Bf°^ ~ [10"-10"] G. We indicate with to the moment at 
which the new internal poloidal component is fully devel- 
oped. 

In order for the poloidal field to affect the spin frequency 
evolution of the star it is necessary that the it diffuses out- 
side of the core. A crucial issue concerns the time-scale Tdiss 
on which the internal poloidal component can diffuse. Un- 
fortunately there are no precise estimates about the value of 
Tdiss in the presence of a layer of superconducting material 
in the inner crust. A detailed description of the diffusion of 
the internal magnetic field is beyond the scope of this pa- 
per. We limit ourselves to note that if Tdiss is shorter than 
a few ten million years the star enters a propeller regime 
too quickly and it is impossible to accelerate the NS to 
the higher observed spin frequencies. On the other hand 
if Tdiss is much larger than lO^'^ yr the expulsion of the 
magnetic field would be phenomenologically irrelevant. It is 
thus necessary that the diffusion of the internal poloidal field 
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Figure 1. (color online) Evolution of the external magnetic 
field (yellow area) assuming a diffusion time-scale t^^^^^ ~ 
[10^—10^''] yr. The blue area indicates the value of the internal 



poloidal component B; 



pol 



[10l°-10"] G which is about two 



orders of magnitude smaller than the toroidal component 5*°*^ 
Here B^°J(to) = 108 G. 



takes place on a time-scale Tdtss ~ [10*-10^°] yr. These esti- 
m ates are compat i ble w ith t he theoretical predictions Rive n 
bv lGeppert et al.l l| 19991 ) and lKonenkov fc Gepper^ l|2000l) . 
We begin our analysis at i = io- We follow the evolution 
of the external magnetic field B^°^ while the internal field 
Bf° diffuses outside of the core. For simplicity we assume 
an exponential growth of the external magnetic field: 



B'::i it) = BLl(io) + Bfr [l - exp(-At/rdi, 



(1) 



where At = t — to and the initial external magnetic field is 
assumed to be _Bg°j(fo) ~ 10* G. The latter is a rather typ- 
ical strength for neutron stars in Low Mass X-ray Binaries 
(LMXBs). 

In Figure [T] we show the evolution of the external mag- 
netic field obtained by using Eq. ([T|. 

In the following we consider the diffusion of the new 
generated internal magnetic field separately in the case in 
which mass accretion is still active and in the case it has 
terminated. 



2.1 Evolutionary scenario without mass accretion 

Here we consider the evolution of a recycled neutron star in 
a scenario in which the mass accretion phase ends within 
a few million years after to. In this case the evolution of 
the spin period of the star is quite simple because it is not 
necessary to take into account the interaction between the 
accretion disk and the magnetic field during the diffusion of 
-Bp";. The total angula r momentum of the star satisfies the 



equation (|BeckeiJl2009l ) 



dJ 



dJ _ It) /7r\ 
m ^ 3" Vc/ 



p3 



(2) 



where c is the speed of light, J = IQ, is the angular momen- 
tum of the star and 3^°^ is obtained from Eq. [T] 



Figure 2. (color online) Evolution diagram of accretion-powered 
X-ray pulsars and of rotation-powered radio pulsars. The data 
points are known rotation-powered pulsars; tho se of pulsars in bi- 
nary systems are encircled. The data are from [Manchester et al.l 
(200S ). The two anomalous stellar objects under study are la- 
beled. The das hed line corresponds to sta rs accreting at the Ed- 
dington limit JLamb &: Boutloukoj 120081) . The shaded area in- 
cludes possible paths (indicated by the red arrows for the case 
without mass accretion) of MSPs after to (see text). Notice that 
the two anomalous stars lie well inside the shaded area. 



In the following we indicate with Po the spin period at 
the time t — to and we treat Po as a free parameter. No- 
tice that the internal magnetic field can develop only if the 
star e nters the r-mode instability region (Cuofano fc Dragol 
I2OI0I ). The upper limit of this region is characterized by 
a period Phigh whose precise value depends on the viscous 
properties of t he star and on the pos sible existence of the 
Ekman layer IjBondarescu et al.l 120071 ). Therefore also the 
parameter Po is bound by the condition Pq $; Phigh ^ 5 ms. 

We stress that Bf° and Pq are not the same for every 
star. They actually depend on the very complicated evolu- 
tion of the internal magnetic field which in turn is related 
to the values of M and B^ °f(to) (see for instance Figure 3 
of ICuofano fc Dragol (|2010l )). 

Results for the evolutionary paths in the P-P plane 
are shown in Figure [21 The shaded strip includes the pos- 
sible trajectories of the stars (the red arrows are included 
as guidelines for the temporal evolution). The shape of this 
strip is regulated mainly by the three parameters Tdiss , Bf° 
and Phigh {B^°f{to) plays a marginal role). The width of the 
strip at the beginning is fixed by Phigh and Piow The lat- 
ter parameter corresponds to the minimal period reachable 
by a recycled MSP. Observationally its value is of the or- 
der o f Plow ~ 1.3 ms (jChakrabartv et al.ll2003l : IChakrabartvl 
I2OO5I 1H Bf° fixes the width of the strip at late times and 
Tdiss regulates the height of the elbow: the smaller Tdiss the 
higher the elbow. 



^ This value can be justified theoretically by considering the 
deformation of the s tar generated by the internal toroidal field 
llCuofano et al.ll2012h . 
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Figure 3. (color online) Panel (a): time evolution of Tco (green 



line) , 



(red line) and period P of the star (blue line). 



The dashed black line indicates the moment at which the spin 
equilibrium phase takes place. Here we assume Pq = 5 ms, 
Bf°^ = lO^O G and M = lO-^Moyr"!. Panel (b): frequencies 
at which the spin equilibrium phase takes place as a function of 
the mass accretion rate. We show two examples, corresponding 
to Pq = 5 ms (gray area) and Pq = 2 ms (cyan area). The width 
of the strips corresponds to the uncertainty on B'L^j (see yellow 
area in Figure [TJ. The red dot refers to the case shown in panel 
(a). 



Remarkably, we can populate the region above the line 
corresponding to the Eddington limit (dashed line). This is 
exactly the region where J1823-3021A and J1824-2452A lie. 



2.2 Evolutionary scenario with mass accretion 

If the mass accretion phase lasts tens or hundreds of million 
years after to, the evolutive scenario is more complicated. 

The plasma flowing onto the neutron star forms an ac- 
cretion disk whose inner edge is given by the coupling radius 
of the magnetosphere (of the same order of magnitud e of the 
Alfve n radius) r^ag ~ B'^''^ R^'^^\My/2GM)-'^''' (iTaurid 
[2013). The interaction between the magnetic field of the 
compact object and the conducting material flowing from 
its companion can provide the nec essary torque Jg = 
M(GMr mag)^^^ to spin up the pulsar jShapiro et al. Ill983h . 
However, under certain conditions, it is possible that this 
interaction leads to a torque reversal. The growth of the ex- 
ternal field -Be°t causes the magnetic boundary rmag to move 
outward relative to the corotation radius r^o = {GM/^^Y'^ 
defined as the distance at which the spin frequency of the 
star is equal to the Keplerian frequency. The star enters a 
propeller phase when rmag > rco: a centrifugal barrier pre- 
vents the material to flow onto the star and a new braking 
torq ue Jm ^ —M fGMr mag)"*^ acts to slow down the pul- 
sar l|Taurisll20lj) . 

In Figure [3^ we show a typical evolutionary path for an 
accreting neutron star if the diffusion of the internal poloidal 
component Bf° is taken into account. The growth of the 
external magnetic field B^°^ leads to the end of the spin-up 
phase. The star evolves into an equilibrium spin phase, it 



alternates between spin-up and spin-down (rmag ~ ^co) and 
finally it appears as a radio pulsar at the end of the mass 
accretion phase. 

Finally, from the analysis presented above we can notice 
that in order to reach high spin frequencies it is necessary a 
mass accretion rate M > 10~^ Moyr"^ (see FigurelSla). 



3 MILLISECOND PULSARS J1823-3021A AND 
J1824-2452A 

In the previous Section we have presented a general scenario 
for the evolution of MSPs when the production and the dif- 
fusion of the internal magnetic field is taken into account. In 
this Section we concentrate on the cases of the two isolated 
pulsars J1823-3021A and J1824-2452A. 

In the most widely discussed model, MSPs have been 
re-accelerated through mass accretion, in particular by 
having spent a fraction of their lives as part of a LMXB. 
The way by which they end up as isolated pulsars is still not 
fully understood. A possibility is that the companions are 
destroyed through abla tion caused by energetic radiatio n 
from the pulsar (Bhat tacharva fc van den Heuvell Il99ll ). 
Another possibility requires mass loss from the companion 
which can take place in very close binary systems due 
to tidal dissipation (JRadhakrishnan fc Shukrd [l98g). A 
totally different scheme to produce an isolated MSP is the 
one in which the pulsar had an interaction with another 
star or binary in the globular cluster, and was spun-up 
in the process (|Lvne et al.lll996l) . In the following we will 
concentrate on the scenario in which MSPs are produced in 
LMXBs. 

The pulsar J1823-3021A, locat ed in the globular cluste r 
NGC 6624, was discovered in 1994 (JBiggs et al.lll990l . ll994l ). 
The Fermi Large Area Telescope ha s recently detecte d the 
7 ray counterpart of its pulsations (IFermi Lat 1120111 ) with 
a luminosity of Lj — 8.4 ± 1.6 x 10'^* ergs s~^, the highest 
observed 7 ray luminosity for any MSP. The pulsation 
period P is 5.44 ms and the observed spinning down rate 
is extremely large Fobs = -1-3.38 x 10~^* s s~^, larger by 
far than the typical value for other MSPs. These properties 
make J1823-3021A an extremely interesting stellar object. 
The total observed 7 ray emission implies that a significant 
fra ction of Poi,s is d ue to the intrinsic spinning down rate 
P (JFermi Lat l l: 



201 ll ). Within the standard magnetic dipole 



model for pulsar emission, such a large P is generated by a 
large surface magnetic field Bo , larger than 10^ G. This fact 
repre sents a puzzle for standard MSP theory (|Fermi Lat I 
I2OIH ): if J1823-3021A was spun-up by mass accretion 
from a companion, due to its large magnetic field Bo, 
even assuming an accretion rate at the Eddington limit 
Msdd = 1Q~^ Mq / yr , the star would have reached a period 
P — 6.6 ms, which is larger than its present period. 

The isolated millisecond pulsar J1824-2452A is located 
in the globular cluster M28 at 5.6 kpc from the Sun 
JLvne et al.lll987l ). The pulsation period is P = 3.054 ms 
and the observed spinning down rate is Poba = -1-1.62 x 10~^* 
s s~^. Indirect evidence exists that this large P observed is 
indicative of the tr u e valu e, e.g. the detecti on of glitches 
(jCognard fc Backed |2004| ). giant pulses JKnight et all 
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[200g), strong X-ray pulses (jRots et al. 1998) and a large 
second period derivative P (Cognard ct al. 1996). As for 
J1823-3021A, the large P implies a surface magnetic field 
-Bo > 10® G and the combination of values (P-P) is difficult 
to explain in the standard MSPs theory. 

Note that these anomalous pulsars can be explained 
in the standard model if the accretion occurred at super- 
Eddington rate. About such a possibility, we note that all the 
accreting compact stars detected up to no w in LMXBs ha ve 
steady mass accretion rates M <C MEdd (|Gallowavi 120081 rl. 
Moreover, taking into account also the interaction between 
the external magnetic field and the accretion disk, a minimal 
value for the mass accretion rate M > 5 x 10~* Mq yr~^ 
is necessary to accelerate the stars up to the present spin 
frequency. For lower values of M the accreting NS enters a 
propeller regime that prevents a further acceleration. 

In our analysis we have presented a scheme in which 
the MSPs can populate the region occupied by J1823-3021A 
and J1824-2452A. It is important to remark that the exis- 
tence of a propeller phase does not invalidate the possibility 
to interpret these two special MSPs in our scheme. Indeed, 
in Figure [3^ we show that the spin frequency of the star 
does not change significantly during the spin equilibrium 
phase: even if a fast spi n-down of the star takes place dur- 
ing the propeller regime (lTaurisll2012l ). the star nevertheless 
will continue to appear as a radio pulsar in the shaded area 
of Figure [2] Concerning the two anomalous stellar objects 
they are no more accreting and if a propeller phase took 
place during accretion it is clearly already over. The exact 
position in the P-P plane before a possible propeller phase 
cannot be determined precisely, but it was located at higher 
frequencies: the existence of this phase would make there- 
fore even more difficult to interpret these two objects within 
the standard pulsar model (because they would lie further 
from the Eddington limit) but it is easy to accommodate in 
our model. 



4 CONCLUSIONS 

We have discussed the evolution of the external magnetic 
field of recycled MSPs due to the diffusion of strong internal 
magnetic fields generated by r-modes during the mass ac- 
cretion phase. Our analysis is based on the previous work of 
ICuof ano & Drago (2010) and Bonanno ct al. (2011) where it 
has been shown that r-modes generate strong internal mag- 
netic fields in accreting compact stars. We have shown that 
the growth of the external magnetic field affects the evo- 
lution of the spin-down rate P of the MSPs and that this 
mechanism can explain the high P of J1823-3021A and of 
J1824-2452A. Our analysis open new evolutionary paths for 
the recycled MSPs in the P-P plane. 



^ LMXBs are often transient sources and may have s hort-term 
active periods of intense accretion; large M presented in lGallowavl 
1120081 ) in parentheses refer to these active periods and not to 
steady accretion. 
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